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摘 要 : 为 探究 根系 水 分 胁迫 响应 函数 对 农田 水 分 动态 及 产量 模拟 的 影响 ,基于 Richards 方 程 和 
PS123 作物 生长 模型 分 别 进行 了 土壤 水 分 动态 和 小 麦 产 量 的 模拟 ,对 比分 析 了 VG(S 型 曲线 )、MP 


(Hy oh a 


Be) KLS(S Æ th A )3 AP KA ME D7 PB KA LE RSE (3 a) 及 潇 河 站 (2 a) 的 试 


验资 料 对 模型 中 的 土壤 水 力 特征 参数 、 水 分 胁迫 响应 函数 参数 进行 反 演 , 确 定 最 优 的 参数 值 ,得 到 


了 土壤 含 


水 率 、 蒸 散 量 及 籽粒 产量 的 模拟 结果 。 结 果 表 明 :(1) VG、MP 及 LS 函数 条 件 下 的 土壤 水 
分 及 产量 模拟 结果 均 较 好 ,回归 估计 标准 误差 (RMSE ) 值 在 0.021~0.036 之 间 ,模拟 值 与 实测 值 间 的 
民 检 验 达 到 了 极 显 著 水 平 ;对 于 霍 泉 站 ,VG .MP LS 函数 条 件 下 ,土壤 含水 率 模 拟 值 与 实测 值 的 平 


均 相 对 误差 分 别 为 6.37% 、8.26% 及 7.18% ,相关 性 系数 值 最 小 为 0.7814, 模 拟 产 量 的 平均 相对 误差 


值 分 别 为 8.73% 、8.40% 及 8.42%; 对 于 潇 河 站 ,土壤 含 


水 率 模拟 值 与 实测 值 的 平均 相对 误差 结果 普 


aA TERI, RAHA 12.47% (2) 土壤 水 分 动态 模拟 ,S 型 曲线 的 水 分 胁 迪 响应 函数 较 凹 凸 型 曲 


线 的 水 分 胁迫 响应 函数 表现 出 了 更 好 的 效果 ,其 


中 使 用 VG 函数 模拟 的 平均 相对 误差 较 小 ,模拟 的 


蒸 散 量 更 接近 于 实测 值 。(3) 产量 模拟 ,3 种 水 分 胁迫 响应 函数 差异 不 明显 。 综 上 ,VG 函数 是 一 种 


精度 较 高 的 根系 水 分 胁迫 响应 函数 , 且 模 型 简洁 方便 ;改进 的 LS 函数 能 够 提高 模拟 精度 ,但 模型 


的 稳定 性 有 待 近 一 步 研究 。 


关 键 词 : 水 分 胁迫 ; 根系 吸水 模式 ; 农田 水 分 模拟 ; 产量 模拟 ; 参数 反 演 


文章 编号 : 


根系 水 分 胁迫 响应 函数 被 定义 为 实际 根系 吸 
水 速率 与 最 大 根系 吸水 速率 之 比 , 是 研究 干旱 胁迫 
对 根系 吸水 影响 的 重要 参数 。 根 系 吸水 是 植物 获 
取水 分 .养分 .进行 光合 作用 的 基础 ,也 是 联系 水 循 
环 和 碳 循环 的 纽带 ，。 在 土壤 水 分 动态 及 产量 的 
模拟 中 ,根系 吸水 模型 是 极为 重要 的 组 成 部 分 ”， 
模拟 效果 直接 影响 土壤 水 分 动态 模拟 及 产量 模拟 
精度 ”。 安 观 条 件 一 下 的 经 验 模型 “常用 于 模拟 根 
系 吸 水 过 程 ,模型 主要 考虑 根系 吸水 同根 系 密度 分 
布 .土壤 剖面 含水 率 AR ERIKA. TEE 
根系 吸水 速率 往往 难于 直接 测试 ,因此 康 绍 忠 等 ” 


收 稿 日 期 : 2021-06-26; fit] ABA: 2021-10-16 


直接 分 析 了 土壤 剖面 充分 湿润 时 的 根系 吸水 速率 
随 深 度 的 变化 规律 ,并 依据 土壤 -植物 的 水 平衡 关 
系 得 出 了 土壤 剖面 均匀 湿润 时 的 根系 吸水 速率 模 
型 ,用 于 计算 最 大 根系 吸水 速率 ,在 实际 应 用 中 较 
为 简单 。 

土壤 剖面 含水 率 对 根系 吸水 有 直接 影响 , 当 土 
培 含 水 率 降 低 时 ,根系 吸水 受到 抑制 ,吸水 速率 小 
于 充分 湿润 时 的 吸水 速率 ”。 根 系 水 分 胁迫 响应 函 
数 用 于 反映 土壤 水 分 胁迫 对 根系 吸水 的 抑制 作用 ， 
使 根系 吸水 模型 可 以 适用 于 干旱 、 半 干旱 条 件 , 朱 
永 华 等 "引入 土壤 湿度 有 关 的 参数 改进 了 根系 吸水 
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模型 ,得 出 了 干旱 荒漠 区 域 的 根系 吸水 模型 。 根 系 
水 分 胁迫 响应 函数 主要 表征 为 土壤 含水 率 或 土壤 
基质 势 的 分 段 函 数 , 有 线性 ”和 非 线 性 "两 类 。 

线性 函数 表明 根系 吸水 速率 随 士 壤 基质 势 或 十 
坏 含 水 率 的 下 降 呈 线性 减 小 趋势 ,其 中 Feddes 4 
所 提出 的 分 段 线性 函数 最 为 典型 。 非 线性 函数 分 
为 S 型 函数 中 和 四 凸 型 函数 52 种 ,二 者 均 表 明 根 
系 吸水 速率 随 着 土壤 基质 势 或 土壤 含水 率 变化 是 
一 个 渐变 的 过 程 ,其 中 分 别 以 Van Genuchten 模 型 "9 
和 康 绍 忠 等 ”建立 的 适用 于 小 麦 根系 吸水 模式 的 指 
数 模型 较为 简便 、 典 型 。 与 线性 函数 相 比 , 非 线 性 
的 土壤 水 分 胁迫 响应 函数 在 大 多 数 情 况 下 能 更 加 
准确 地 表征 土壤 水 分 对 根系 吸水 的 影响 下。 四 是 
型 响应 函数 在 分 段 连 接 处 会 出 现 变 化 速率 突变 ,为 
避免 该 现象 本 研究 提出 了 以 土壤 含水 率 为 变量 的 S 
型 水 分 胁迫 响应 函数 模型 ,与 男 外 2 种 典型 的 根系 
水 分 胁迫 响应 函数 进行 对 比分 析 ,该 对 比 在 作物 水 
模拟 过 程 中 完成 。 

作物 水 模型 主要 用 于 对 作物 生长 发 育 动 态 、 产 
量 进行 模拟 ,是 定量 化 评价 灌水 施肥 等 田间 管理 措 
施 的 重要 手段 " ,其 模拟 对 象 主要 为 土壤 水 .土壤 
氮 素 和 作物 生长 。 土 壤 水 的 模拟 主要 是 土壤 水 分 
运动 与 蒸 散 过 程 , 不 同 的 模型 采用 方法 不 同 ,如 在 
水 分 模拟 中 ,WOFOST"” DSSAT®” PS12322 等 模型 
采用 计算 较为 简单 的 水 平衡 法 ,而 HYDRUS-1D*、 
DAISY 等 模型 则 采用 了 以 Richards"" 方 程 为 基础 
的 土壤 水 动力 学 模型 ,更 适用 于 边界 条 件 复 杂 的 环 
境 且 计算 精度 较 高 ” 。 作 物 蒸 散 模拟 中 大 部 分 模 
型 采用 了 Penman-Monteith 方程 计算 参考 作物 蒸 散 
量 , 通 过 作物 系数 计算 潜在 蒸 散 量 ,最 后 使 用 叶 面 
积 指数 将 潜在 蔡 散 量 分 为 潜在 芋 腾 量 和 潜在 蒸发 
量 , 实 际 蒸发 量 可 作为 Richards 方 程 的 上 边界 条 件 ， 
潜在 蒸腾 量 输入 作物 根系 吸水 模型 作为 实际 根系 
吸水 的 驱动 。 作 物 生长 模拟 过 程 是 从 光合 有 效 辐 
射出 发 ,考虑 光 、 温 、 水 和 养分 等 胁迫 ,描述 光合 作 
用 及 干 物质 积累 分配 和 转化 过 程 。 目 前 , 既 有 较 
通用 的 WOFOST , DSSAT DAISY 等 模型 ,也 有 针对 
不 同 作物 的 特定 模型 ,如 用 于 模拟 谷类 作物 的 CE- 
RES27 模 型 。 以 荷兰 De Wit 教 授 为 首 的 研究 小 组 ， 
先后 建立 和 完善 了 在 不 同 生 产 力 水 平 下 的 作物 生 
长 模拟 模型 ,以 PS123 模 型 为 代表 ,更 多 地 强调 作物 
的 共性 ,与 作物 的 生理 生态 机 理 结合 较 紧 , 在 我 国 
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的 适应 性 方面 已 有 良好 的 工作 基础 “, 如 毛 振 强 等 ” 
以 河北 曲 周 县 30 a 的 气象 数据 为 依据 ,应 用 PS123 
模型 分 析 了 该 地 区 夏 玉 米 的 需 水量 和 优化 灌溉 方 
案 , 对 指导 当地 玉米 生产 具有 重要 意义 。 

PS123 模 型 是 联系 水 分 与 作物 生长 关系 的 主要 
环节 ,其 对 于 水 分 动态 模拟 采用 简化 的 水 平衡 法 ， 
对 根系 吸水 只 是 依据 机 理 进行 了 较为 简单 的 线性 
化 处 理 , 需 要 进一步 完善 ”。 鉴 于 此 ,本 研究 选用 
PS123 作物 生长 模型 与 Richards 方 程 建立 了 作物 产 
量 与 水 分 关系 ,并 采用 机 理性 较 强 的 根系 吸水 模 
型 。 采 用 规划 求解 方法 对 模型 参数 进行 率 定 ,对 不 
同形 式 水 分 胁迫 响应 函数 的 表征 效果 进行 了 评价 。 


1 研究 区 概况 


1.1 试验 站 概况 

以 冬小麦 作为 研究 对 象 ,以 山西 霍 泉 、 潇 河 两 
处 试验 站 的 试验 资料 为 依据 。 霍 凡 站 位 于 霍 录 灌 
区 中 上 游 的 洪 洞 县 广 胜 寺 镇 东安 村 (111°46'21”E， 
36"17'35"N ) ,海拔 高 程 529 m ,该 地 所 属 气候 为 暧 温 
带 大 陆 性 气候 ,土壤 质地 为 轻 壤 土 ,0~100 cm 土壤 
容重 1.46 gcm ,田间 持 水 率 为 0.3592 cm cm >。 
清河 试验 站 位 于 晋 中 市 榆 次 区 杨 村 (112*41'31" 开 ， 
37°39'06"N ) ,海拔 高 程 782.64 mm, 土质 为 重 坏 十 ,0~ 
100 cm 土壤 容重 为 1.42 g- cm 3, 田间 持 水 率 为 
0.4044 cm? cm 。 

模拟 中 取 用 了 霍 泉 站 2017 一 2019 年 与 潇 河 站 
2004 一 2005 年 共 5 个 站 年 的 冬小麦 试验 资料 (小 麦 
生育 期 跨越 冬季 ,本 文中 年 为 收获 年 ) ,主要 包括 各 
年 的 气象 数据 .冬小麦 生长 发 育 期 的 干 物质 测试 数 
据 及 土壤 含水 率 实测 数据 。 其 中 ,土壤 含水 率 测试 
方法 , 霍 泉 站 为 烘 干 法 , 潇 河 站 为 中 子 仪 法 。 霍 泉 
站 为 大 田 试验 ,小 区 面积 为 66.7 m, FATED RR T 
无 底 测 坑 , 测 坑 面积 为 3.33 mx2 mo 
1.2 试验 处 理 选 定 

霍 录 站 与 淆 河 站 均 为 冬小麦 水 肥 耦 合 试验 。 
其 中 , 霍 泉 站 每 年 设置 6 个 处 理 , 潇 河 站 每 年 设置 10 
个 处 理 。 本 研究 仅 在 每 年 的 试验 处 理 中 选取 了 高 
水 .中 水 和 零 水 处 理 用 于 模型 参数 的 率 定 ,详细 处 
理 信 息 见 表 1。 由 表 可 见 ,不 同年 份 . 处 理 间 施肥 数 
量 有 所 不 同 ,因此 研究 中 土壤 含水 率 、 蔡 散 量 和 产 
量 的 变化 为 水 肥 双 重 影响 下 的 结果 。 


202204.00136v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


FHEWY 


表 1 冬小麦 水 肥 耦 合 试 验 处 理 基本 信息 


Tab. 1 Basic information of winter wheat water and fertilizer coupling test treatment 


灌水 定额 灌水 时 间 ( 月 -日 ) 底肥 追肥 
站 名 年 份 。 处 理 WK = = T 7 oe 
ea 2017 高 水 75 11-22 03-22 05-08 05-22 600 450 
中 水 75 11-22 03-22 05-22 = 600 300 
零 水 - - - - - 600 - 
2018 高 水 75 03-15 04-18 05-12 05-31 900 563 
中 水 75 03-15 04-18 05-12 - 600 375 
AIK - z - - 7 Š = 
2019 高 水 75 03-26 04-22 05-22 一 1050 510 
中 水 75 03-26 05-08 = 一 750 345 
零 水 = z - = = - = 
潇 河 站 2004 高 水 75 12-11 04-19 05-07 05-29 750 225 
中 水 45 04-19 一 一 一 750 225 
零 水 = = - - - 750 - 
2005 高 水 75 12-01 04-02 05-11 05-31 1125 338 
中 水 45 04-20 = 一 一 750 225 
Zk = = = - - 750 - 
dK(0.) 
2 ”研究 方法 TN (6) 
2.1 土壤 水 热 耦合 变化 动态 模拟 D(0)=K(O) Sr) ie 
2.1.1 土壤 水 热 耦 合 动态 模拟 模型 ”本 研究 在 冬 小 jp Re iek 
麦 生 育 期 水 分 模拟 的 过 程 中 考虑 了 温度 变化 ,采用 a(h)=4 [ s(a ` (8) 
包含 根系 吸水 项 的 Richards 方程 进行 土壤 水 分 动态 6. h>h, 
模拟 ,以 1d 作 为 模拟 时 段 ,对 距 地 表 2 m 内 的 土壤 下 a)? 
含水 率 及 热流 进行 分 层 模拟 (20 em 一 层 ) ,依据 kep- K3) UES j 0<0, (9) 
STVF 理论 模式 ”分 析 土 壤 温度 变化 对 土壤 水 分 运 a a 
K, 0>0. 


动 特性 参数 的 影响 。 采 用 尚 松 浩 等 ”给 出 的 水 热 
耦合 方程 模拟 热流 变化 ,该 方程 综合 考虑 了 土壤 水 
分 迁移 .热传导 水 分 相 变 等 因素 ,相对 于 传统 水 热 
耦合 方程 ,水 . 热 方程 间 的 耦合 性 大 为 减少 ,可 以 不 
考虑 相 变 作用 ,只 在 计算 时 段 末 根 据 冻 结 曲线 进行 
冰 和 未 冻 水 含量 的 修正 ,有 利于 提高 方程 求解 的 挝 
代 计 算 。 


30 _ a0 dg] 9K(8) _ 
0 -30 ‘D(6) £] = Sle) (1) 
T al OT eat 
C. ðt = 二 .| U. dz (2) 
do 
C=Lip, IT (3) 
C.=C,+C, (4) 


A,=A+DEO)C, (5) 


式 中 :6 为 土壤 体积 含水 率 (cm + em™*) ;1 为 时 间 
(min) ;D(9) 为 土壤 非 饱 和 扩散 率 (cm + min”) ,根据 
式 (7) 计 算 ,其 中 怨 为 比 水 容量 ,9 与 4 的 关系 由 土 


壤 水 分 特征 曲线 确定 ( 式 8) 89s KON BEAK 
和 导 水 率 (cm:min :) ,采用 目前 普遍 使 用 的 式 (9) 计 
算 …225;z 为 深度 (cm) ;SGz 思 为 根系 吸水 速率 (min 7 ); 
7 为 土壤 温度 (%C);L 为 融化 潜 热 ,研究 中 取 335 J+ es; 
0, 为 相应 土壤 负 温 度 条 件 下 可 能 的 最 大 未 冻 水 含量 
(em*cm”) ;0, 为 土壤 容积 比 热 (J.g'…%C"') ,C=1.92 
(1-0.,)+4.189 ;A 为 土壤 体积 热 导 率 (Wcm °C"), 
A=C.D,, DJ ERAP CK , D=0.002+0.01150+ 
0.01030%; 0 为 土壤 未 冻 水 体积 含水 率 (cm cm- ) ， 
未 冻 水 含量 可 以 利用 脉冲 核磁 共振 法 . 量 热 法 .时 
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域 反射 法 等 方法 测定 " ,其 测定 过 程 较为 复杂 ,本 
文 直接 引用 了 尚 松 浩 给 出 的 冻结 特征 曲线 ;90. 为 
饱和 土壤 含水 率 (cm: cem’); C C A UAIR 
的 相 变 热 容 量 、 等 效 体积 热 容量 、 等 效 热 导 率 及 等 
效 对 流速 度 , 由 0,、9,、A、D(9)、K(9) 通 过 公式 计算 得 
到 。0. 为 残余 土壤 含水 率 (cm:cm’); 有 为 土壤 基质 
势 (cm);an 为 土壤 水 分 特征 曲线 形状 参数 ,其 值 通 
过 参数 反 演 确定 ;m=1-1/n;h. 为 经 验 参 数 (cm) ,本 
人 研究 取 值 为 0;K. 为 土壤 饱和 导 水 率 (cm:min ') ,其 
值 通 过 参数 反 演 确定 ;6, 为 模型 曲线 外 插 参 数 ,数值 
比 9. 稍 大 。 

上 述 水 热 耦 合 方程 与 初始 条 件 .边界 条 件 一 起 
组 成 完整 的 水 热 耦 合 动态 模拟 模型 。 初 始 条 件 A 
界 条 件 采用 王 仰 仁 呈 论 文中 给 出 的 计算 方法 ,初始 
条 件 中 土壤 含水 率 为 播种 日 土壤 剖面 含水 率 ,气温 
根据 自动 气象 站 测 得 ,地温 采 用 地 温 温度 计 测 得 ; 
上 边界 条 件 中 地 表 通 量 等 于 棵 间 土 壤 蔡 发 量 '”, 地 
表 温 度 根据 空气 温度 与 地 表 温 度 的 经 验 关系 计算 ; 
下 边界 条 件 中 的 温度 及 土壤 含水 率 因 在 小 麦 生育 
期 内 变化 较 小 ,因而 采用 第 一 类 边界 条 件 ”。 利 用 
差分 方法 求解 ,可 得 分 层 土壤 含水 率 及 温度 随时 间 
的 变化 过 程 。 
2.1.2 根系 级 水 模型 ”采用 潜在 的 根系 吸水 速率 与 
水 分 胁迫 响应 的 乘积 计算 根系 的 实际 吸水 速率 ”1。 

S(z,l)=S0(z,)B (10) 

式 中 .Ss(z,) 为 根系 吸水 速率 (min') ;So(z 思 为 土壤 充分 
湿润 条 件 下 的 根系 吸水 速率 (min ') , RH Novák” 2A 
出 公式 计算 , 见 式 (11);pB 为 水 分 胁迫 响应 函数 ,是 
以 土壤 基质 势 或 土壤 含水 率 为 自 变 量 的 函数 。 


S (21) = Tp() (11) 


zh = e] 

式 中 ;Tp(?) 为 作物 潜在 蒸腾 速率 (cm:min') ,是 作物 
潜在 蒸 散 量 的 一 部 分 ,可 根据 经 验 公式 计算 得 到 ; 
z 为 深度 (cm ) ;zs 为 根系 层 深度 (cm) ;6 为 经 验 系数 ,与 
作物 种 类 有 关 , 本 文 参照 由 懋 正 ” 的 研究 取 6=3.8。 

本 研究 对 3 种 形式 的 水 分 胁迫 响应 函数 进行 了 
对 比分 析 。 

(1) 以 土壤 基质 势 为 因 变 量 的 Van Genuchten 
函数 中 (以 下 称 为 VG)。 

B= 一 (12) 
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式 中 :BO 为 以 土壤 基质 势 为 自 变 量 的 水 分 胁迫 响 
应 函数 , 呈 S 形 状 变化 ;h@) 为 地 面 以 下 深度 z 处 的 土 
壤 基 质 势 (cm) ;7(z) 为 土壤 溶质 势 , 本 研究 未 考虑 溶 
质 影 响 , 因 此 取 mw(z)=0;P 为 经 验 系 数 ,通过 参数 反 演 
确定 ;js 为 作物 蒸腾 量 减 小 到 最 大 可 能 蒸腾 量 50% 
时 所 对 应 的 土壤 水 势 (em) ,其 值 通过 参数 反 演 确定 。 

(2) 以 土壤 含水 率 为 因 变 量 的 需 函 数 (以 下 
称 为 MP)。 


1 0>0, 
B0) = (AW) 0,,<0<0, (13) 
0 0<0,, 
Awe Ow (14) 
6-0 


式 中 :B(9) 为 以 土壤 含水 率 为 自 变 量 的 水 分 胁迫 响 
应 函数 ,时 四 ( 凸 ) 形 状 变 化 ;AW 为 有 效 土 壤 含 水 率 
的 相对 值 ;0. 为 z 深 上 度 的 土壤 含水 率 (cm:cm”);0 为 
临界 土壤 含水 率 (cm:cm’);0,, 为 凋 茸 土壤 含水 率 
(em-em™) ;A HAIMA R 5 0.0. AA HBB 
确定 。 

(3) 以 土壤 含水 率 为 因 变 量 的 改进 型 分 段 函 
数 , 式 (13) 在 分 段 连 接点 处 会 出 现 速率 突变 ,不 符 
合作 物 生 长 的 连续 渐变 过 程 ,鉴于 此 引入 参数 B 、C 
消除 速率 突变 ,形成 $ 型 曲线 ,由 此 给 出 了 LS ,计算 
公式 如 下 : 


1 0>0 
ol 
“o= 1 
BO)= b+] 0 S00, (15) 
0 0<0 


式 中 : BO) 为 改进 线 型 的 水 分 胁迫 响应 函数 ;B.C 为 
经 验 系 数 ,通过 参数 反 演 确定 。 
2.2 冬小麦 产量 模拟 

采用 PS123 作物 生长 模型 进行 冬小麦 产量 模 
拟 ,模拟 过 程 中 除了 考虑 有 效 太阳 辐射 和 温度 外 ， 
同时 考虑 有 效 水 分 和 有 效 养分 胁迫 的 影响 。 模 拟 
产量 利用 大 田 作物 总 的 同化 潜力 表示 ” ,计算 过 程 
如 下 : 

Fgass, = Fgc, x 30/44 x FW, XFN,XCVF (16) 

式 中 :Fgass; 为 大 田 作物 总 同化 率 (kg'hm”.d);Fgc， 
为 密闭 参照 作物 总 CO: 吸 收 率 (kg'hm”d ) ,根据 日 
长 温度 . 叶 面 积 指数 .消光 系数 . 光 能 利用 效率 、 冠 
层 项 部 截获 的 光合 有 效 辐 射 计算 ;30/44 为 CHO 
与 C0,; 的 分 子 量 之 比 ;FW 为 水 分 胁迫 系数 ;FN 为 养 
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分 胁迫 系数 ,分 别 采用 式 (17) ”和 式 (18) 计 算 ;CVF 
为 光合 产物 转化 效率 ,通过 参数 反 演 得 到 ;i 为 生长 
天 数 。 


A, 
FW,= Ey a (17) 
i 7 P; 


Sr, > Sr; 
ml Bn 
FN, = i i (18) 


式 中 :ETa 为 水 分 胁迫 下 的 作物 实际 蒸 散 量 (mm.d); 
ETp, 为 作物 潜在 蒸 散 量 (mm.d7) ;Sr 为 根系 吸 氮 量 
(kg"hm*…d') ,根据 实际 根系 吸水 量 与 土壤 水 溶液 
氮 素 浓度 的 乘积 计算 中 ;Sm 为 植株 最 低 限 度 含 氛 量 
(kg*hm ”dd ) ;A 分别 为 水 分 亏 缺 敏感 指数 和 养 
分 亏 缺 敏感 指数 ;hh 三 为 经 验 系数 , 均 通 过 参数 
反 演 得 到 。 
2.3 参数 反 演 

本 研究 利用 Excel 规划 求解 工具 中 的 演化 算法 
反 演 确定 待 求 参数 。 对 霍 果 站 ,将 3a 共 9 个 处 理 数 
据 放 在 一 起 进行 参数 反 演 ;对 潇 河 站 ,将 2a 共 6 个 
处 理 数据 放 在 一 起 进行 参数 反 演 。 反 演 参 数 主 要 
包括 土壤 水 分 动态 模拟 中 的 参数 和 小 麦 产 量 模拟 
中 的 参数 。 

(1) 土壤 水 分 动态 模拟 参数 反 演 ,以 实测 土壤 
含水 率 及 模拟 土壤 含水 率 误差 平方 和 最 小 为 目标 
国 数 , 见 式 (19 ): 


m Pj E 
min > SW; sens S00 | (19) 
j5 r= 


i 


式 中 :SW 为 第 j 个 处 理 所 有 模拟 土壤 含水 率 与 实测 
土壤 含水 率 的 误差 平方 和 ;0; 为 实测 土壤 含水 率 
(cm cm ); 05 为 模拟 土壤 含水 率 (cm cem”); p Hj 
处 理 中 实测 土壤 含水 率 的 个 数 ,i=1, 2,… ,pj;j=1， 
2,…,m, 其 中 以 为 规划 求解 所 用 到 处 理 个 数 , 霍 呆 
站 m=9, 潇 河 站 m=6。 

待定 参数 包括 土壤 水 分 特性 参数 0.、0.、a、n K 
以 及 水 分 胁迫 参数 hso、P、6,,、9,、.4、B、C。 演 化 时 , 根 
据 已 有 的 经 验 参 数 ”… ,参照 试验 站 的 土壤 质地 确定 
3 组 参数 初始 值 ,依次 进行 演化 计算 ,可 获得 3 组 竺 


定 参 数 , 选 定 其 中 土壤 含水 率 拟 合 误差 最 小 的 一 组 
参数 值 作为 初始 值 , 再 进行 演化 计算 , 当 参 数 不 再 
变化 时 停止 演化 计算 。 

(2) 冬小麦 生长 动态 及 产量 模拟 参数 反 演 ,以 
E . 叶 籽粒 模拟 值 与 实测 值 误差 平方 和 最 小 为 目 
标 函 数 , 见 式 (20): 


min >SY 
i=l 


SY, = ŞE, 一 YE’, + Su, 一 Yy + So, -YL' i 
(20) 
式 中 ;SY 为 第 ;个 处 理 中 实测 干 物 质 重 与 模拟 干 物 
质 重 的 误差 平方 和 ;YE YJ YLA AS A, 
籽粒 干 重 (kg.hm2);YE' YJ’ YL’ 分别 为 模拟 叶 、 
25 FPFE (kg hm”) ;egj.4 分 别 为 处理 中 叶 、 
茎 籽粒 实测 干 重 样本 的 个 数 ,o=1, 2,…, e, g=1， 
‘ “+, gvl, 2, tee, lo 
RSRK NI RAAR C kA) FE 
Ay Wh il Be BS BE (hs A2) RIGA PA W ER A eK 
(CVF)。 为 了 获得 更 为 准确 的 产量 模拟 结果 ,在 式 
(20) 的 基础 以 产量 模拟 值 与 实测 值 误差 平方 和 最 
小 为 目标 ( 式 21) ,进一步 进行 参数 反 演 分 析 。 


YL= miny> (YL; - Yuy (21) 
j=1 


X 


N 


式 中 : YL 为 处 理 j 收 获 时 实测 产量 (kg: hm); YL", 
为 /处 理 模拟 产量 (kg'hm”)。 
2.4 拟 合 精度 评价 

采用 回归 估计 标准 误差 (RMSE) 相对 误差 
(RE) 相关 性 系数 (r) 5 对 土壤 含水 率 的 模拟 结 
进行 评价 ,r 值 越 大 ,表明 模型 拟 合 效果 越 好 ， 
RMSE .RE 值 越 小 ,表明 模拟 值 的 精度 越 高 。 采 用 
相关 性 系数 对 地 上 干 物 质 的 拟 合 情况 进行 评价 ， 
收获 时 产量 数据 较 少 ,采用 相对 误差 对 模拟 产量 结 
果 进 行 评价 , 并 利用 检验 对 不 同 水 分 胁迫 响应 函 
数 模 拟 的 产量 进行 差异 性 检验 。 


3 结果 与 分 析 


3.1 水 分 胁迫 响应 函数 对 0~60 cm 土壤 水 分 动态 模 
拟 的 影响 
通过 式 (19) 进 行 参数 反 演 优化 , 取 最 优 演化 结 
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果 , 见 表 2、 表 3。 由 表 可 见 ,不 同 的 水 分 胁迫 响应 函 
数 条 件 下 ,利用 试验 资料 反 演 得 到 的 土壤 水 力 特 征 
参数 .水 分 胁迫 响应 函数 参数 ,没有 明显 的 差异 ,日 
均 在 合理 范围 内 ””, 如 4. 均 大 于 田间 持 水 率 , 晶 不 超 
过 0.5, 表 明了 反 演 参数 的 合理 性 。 

利用 反 演 的 参数 ,可 以 得 到 不 同 土壤 含水 率 模 
拟 结果 ,对 模拟 值 和 实测 值 进 行 RMSE .RE 和 7 的 对 
比分 析 ,结果 见 表 4 及 图 1。 由 图 表 可 见 ,不 同 水 分 
胁迫 响应 函数 条 件 下 模拟 结果 的 显著 性 检验 均 达 
到 极 显著 水 平 ,相关 性 系数 + 最 小 为 0.7196,RMSE 
在 0.0217~0.0363 之 间 ,表明 土壤 含水 率 模 拟 结果 与 
实测 值 具 有 非常 好 的 一 致 性 。 从 霍 泉 站 土壤 含水 


表 2 土壤 水 力 特征 参数 反 演 结果 
Tab.2 Inversion results of soil hydraulic 


characteristic parameters 


站 名 ”模型 g/cmi:cm” 8/cm tem? aem” n KJem+d' 
霍 果 站 VG 0.4468 0.0066 0.0076 1.3769 9.3486 
MP 0.4460 0.0033 0.0245 1.6078 5.5227 
LS 0.4833 0.0095 0.0073 1.8180 2.7249 
潇 河 站 VG 0.4900 0.0036 0.0128 1.4325 7.5811 
MP 0.4502 0.0054 0.0279 1.4296 5.9857 
LS 0.4618 0.0040 0.0303 1.5368 5.1511 


注 :VG .MP LS 分 别 为 以 土壤 基质 势 为 因 变 量 的 Van Genuchten PKI 
数 以 土壤 含水 率 为 因 变 量 的 震 函 数 .以 土壤 含水 率 为 因 变 量 的 改 
进 型 分 段 函 数 ;4. 为 饱和 土壤 含水 率 ;0 为 残余 土壤 含水 率 ;a.n 为 
土壤 水 分 特征 曲线 形状 参数 ;K. 为 土壤 饱和 导 水 率 。 下 同 。 


表 3 3 种 水 分 胁迫 响应 函数 参数 反 演 结果 


Tab. 3 Three kinds of moisture stress response function parameter inversion results 


; VG MP LS 
站 名 
P hso/em Bilem cm lem em” A bcm cm bycm cm B C 
霍 泉 站 1.98 1438.21 0.3586 0.1495 0.5931 0.3709 0.1680 0.5457 2.9801 
漠河 站 1.71 1059.34 0.2574 0.1781 0.5954 0.3781 0.1513 0.5427 2.6097 


TE : ho WIEDZE S ED] BM EAC AY PEZES E 50% 时 所 对 应 的 土壤 水 势 ;0 为 临界 土壤 会 水 率 ;b 为 凋 蓉 土壤 含水 率 ;P.4 .B、C 为 经 验 系数 。 


表 4 不 同 水 分 胁迫 响应 函数 模拟 效果 的 评价 
Tab. 4 Evaluation of simulation effects of different 


moisture stress response functions 


站 名 “模型 RMSE RE 范围 /% 平均 RE/% r 样本 数 
ERR VG 0.0217 0.05~23.82 6.37 0.8557 191 
MP 0.0273 0.01~22.02 8.26 0.8150 191 
LS 0.0223 0.04~27.27 7.18 0.7814 191 
漠河 站 VG 0.0363 0.13~38.46 12.47 0.8280 122 
MP 0.0279 0.08~30.98 8.59 0.7820 122 
LS 0.0314 0.09~40.62 10.33 0.7196 122 


注 :RMSE 为 回归 估计 标准 误差 ;RE 为 相对 误差 ;r 为 相关 性 系数 ， 
检验 结果 均 达 到 极 显 车 水 平 。 下 同 。 


率 模拟 结果 看 ,VG 最 好 ,其 次 为 LS,MP 相对 较 差 ， 
表明 S 型 胁迫 响应 函数 的 模拟 结果 优 于 止 凸 型 琐 
数 ; 对 于 S 型 曲线 ,基于 土壤 基质 势 的 VG 要 好 于 
LS。 从 潇 河 站 模拟 结果 看 ,土壤 含水 率 相 对 误差 大 
TÆR ,平均 RE 最 大 为 12.47%; 对 比 3 种 胁迫 响 
应 函数 ,MP 计算 的 RE 数值 最 小 ,另外 2 种 函数 对 应 
的 相对 误差 均值 较 大 ; 潇 河 站 土壤 含水 率 相 对 误差 
变化 范围 较 大 ,主要 原因 可 能 是 中 子 仪 测定 土壤 含 
水 率 过 程 中 存在 率 定 误 差 ,成 钰 中 的 研究 中 也 出 现 


了 类 似 情况 。 
3.2 水 分 胁迫 响应 函数 对 蒸腾 量 和 蒸 散 量 模拟 的 
影响 


水 分 胁迫 响应 函数 对 蒸腾 量 (Ta) 及 蒸 散 量 
(ETa) 的 模拟 有 明显 的 影响 ( 表 5)。 由 表 5 可见 ,使 
用 VG 模 拟 的 根系 吸水 量 及 ETa 明显 高 于 另外 2 种 
水 分 胁迫 响应 函数 条 件 下 的 模拟 值 , 旦 更 接近 实测 
值 ,2 个 试验 站 非常 一 致 ,不 同 处 理 也 表现 出 了 相同 
结果 。 这 里 ETa 的 实测 值 是 利用 水 量 平 衡 法 计算 求 
得 ,计算 过 程 中 没有 考虑 水 分 的 深层 渗 漏 及 深层 土 
壤 水 的 补给 。 对 于 高 水 .中 水 的 处 理 ,灌水 后 根系 
层 含水 率 较 高 ,水 分 向 下 运 移 的 数量 相对 较 多 ， 
而 灌水 处 理 实测 的 ETa 普 遍 大 于 模拟 值 ;不 灌水 情 
况 下 ,根系 层 含水 率 较 低 ,对 水 分 的 吸 持 力 较 大 ,更 
容易 发 生 深层 土壤 水 的 补给 ,因而 霍 泉 站 零 水 处 理 
模拟 蒸 散 量 普 遍 大 于 实测 值 。 

XIF MPLS 2 种 水 分 胁迫 响应 函数 , 霍 果 站 采 
用 MP 模拟 的 Ta 及 ETa 较 低 ; 而 浪 河 站 采用 MP 模拟 
结果 明显 大 于 LS 模拟 值 。 说 明 MP、LS 受 试验 环境 
变化 影响 较 大 ,存在 不 稳定 性 ,相对 而 言 VG 使 用 效 
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因 变 量 的 Van Genuchten 函数 .以 土壤 含水 率 为 
因 变 量 的 改进 型 分 段 函 数 。 下 同 。 
基于 3 种 水 分 胁迫 响应 函数 的 土壤 含水 率 模拟 值 与 实测 值 对 比 


0.10 


实测 值 /cm?:cm3 


wW 


ZE Er HITE PRA, 


Fig. 1 Comparison of simulated and measured water content values of three moisture stress response functions 


#5 不 同 水 分 胁迫 响应 函数 下 模拟 的 蒸腾 量 及 蒸 散 量 


Tab. 5 Simulated transpiration and evapotranspiration under different moisture stress response functions /mm 
站 名 年 从 模型 
Ta ETa Ta ETa Ta ETa 
霍 果 站 2017 VG 264.4 401.8 224.2 350.3 160.6 263.5 
MP 201.6 312.9 170.2 275.5 111.6 186.7 
LS 244.4 356.7 194.3 297.7 112.2 193.4 
实测 - 502.0 一 400.2 一 261.3 
2018 VG 248.2 383.3 241.7 359.7 150.1 246.7 
MP 186.2 283.3 184.0 266.8 101.0 163.0 
LS 226.1 331.9 221.5 310.6 104.8 176.6 
实测 - 446.4 一 378.7 - 232.7 
2019 VG 245.9 371.7 216.7 323.5 129.9 224.1 
MP 182.6 271.5 157.4 233.5 85.7 150.7 
LS 220.5 315.1 182.0 262.9 82.8 156.0 
实测 一 404.9 一 374.8 一 235.1 
潇 河 站 2004 VG 303.1 451.5 187.5 283.0 130.9 204.0 
MP 270.2 391.8 146.4 224.2 92.0 157.2 
LS 247.8 372.6 128.2 208.9 78.8 149.0 
实测 = 529.6 一 297.3 一 235.6 
2005 VG 337.1 461.1 187.3 234.7 162.1 206.0 
MP 299.3 397.1 141.3 187.8 114.6 157.9 
LS 276.0 380.1 126.1 180.6 96.0 145.9 
实测 一 515.0 一 237.9 一 213.3 


注 :Ta 为 蒸腾 量 ;ETa 为 蒸 散 量 。 
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王 铁 英 等 : 根系 水 分 胁迫 响应 函数 对 土壤 水 及 作物 生长 动态 和 产量 模拟 影响 的 研究 


果 好 。 
3.3 水 分 胁迫 响应 函数 对 产量 模拟 的 影响 

采用 PS123 作物 生长 模型 模拟 作物 产量 ,以 霍 
泉 站 实测 的 地 上 干 物质 重 数据 反 演 了 作物 生长 动 
态 模 拟 参 数 , 并 给 出 了 地 上 干 物质 模拟 值 与 实测 值 
的 拟 合 程度 检验 结果 ,给 出 了 最 终 产 量 模 拟 结果 的 
相对 误差 ( 表 6)。 由 表 6 可 见 ,地 上 干 物质 重 模 拟 值 
与 实测 值 的 相关 性 系数 均 在 0.91 以 上 ,r 检 验 达到 
了 极 显 著 水 平 ,表明 模拟 值 与 实测 值 具有 很 好 的 一 
致 性 ,3 种 水 分 胁迫 响应 函数 都 可 以 获得 高 精度 的 
模拟 结 

由 表 6,VG、MP 及 LS 的 平均 RE 分 别 为 8.73%、 
8.40% 及 8.42% ,产量 模拟 结果 较为 接近 ,图 2 给 出 
了 模拟 产量 与 实测 产量 对 比 情 况 。 由 图 可 见 ,3 种 
水 分 胁迫 响应 函数 模拟 值 与 实测 值 吻合 程度 较 好 ， 
确定 性 系数 在 0.73 以 上 ,模拟 点 贴近 1:1 直线 ,表明 
产量 模拟 结果 较 好 。 对 27 个 模拟 产量 样本 按 不 同 
的 水 分 胁迫 响应 函数 进行 单 因素 的 下 检验 ,计算 下 
IEX 0.0032, /h F Fooi (5.61) ,表明 根系 水 分 胁迫 响 
应 函数 确定 性 对 最 终 产量 模拟 没有 明显 的 影响 。 


4 讨论 


水 分 是 制约 干旱 区 作物 生长 的 主要 因素 ,根系 
水 分 胁迫 啊 应 函数 是 根系 吸水 模型 中 的 重要 组 成 


部 分 。 通 过 评 佑 不 同形 式 水 分 胁迫 响应 函数 的 表 
征 效果 ,有 助 于 完善 对 根系 吸水 过 程 的 处 理 。 在 作 
物产 量 水 分 关系 中 ,根系 水 分 胁迫 响应 函数 通过 以 
下 过 程 影响 土壤 水 动态 、 作 物 生 长 动态 和 产量 的 
模拟 。 

利用 参考 作物 蒸 散 量 、 作 物 系 数 、 叶 面积 指数 
可 计算 潜在 蒸腾 量 潜在 蒸发 量 。 其 中 ,潜在 燕 发 
量 通 过 上 边界 条 件 影响 土壤 水 分 运动 ,潜在 蒸腾 量 
通过 根系 吸水 影响 实际 蒸腾 量 ,同时 根系 吸水 量 也 
用 于 计算 根系 吸 氮 量 。 水 分 胁迫 响应 函数 形式 不 
同时 ,模拟 的 根系 吸水 效果 不 同 , 如 吴 训 等 的 研 
究 表明 , 非 线性 胁迫 响应 函数 相对 于 线性 胁迫 响应 
函数 有 更 好 的 表征 效果 ,本 文 进一步 对 比 了 非 线 性 
水 分 胁迫 响应 函数 间 表 征 效果 ,更 好 地 反映 了 作物 
生长 连续 渐变 过 程 。 土 壤 基质 势 不 受 土壤 分 层 影 
啊 , 随 深度 变化 是 连续 的 ,本文 研 究 结果 表明 VCG 模 
拟 水 分 变化 优 于 男 外 2 种 函数 。 

作物 吸收 的 水 分 .养分 对 生长 发 育 及 生物 量 积 
累 具 有 显著 影响 ,对 此 PS123 模型 将 不 同 影响 因子 
分 为 了 不 同 的 生产 潜力 水 平 ” ,可 采用 实际 蒸腾 量 
与 潜在 蒸腾 量 的 比值 作为 作物 生长 水 分 胁迫 因子 ， 
采用 根系 吸 氮 量 与 植株 最 低 限 度 含 氮 量 的 比值 作 
为 作物 生长 养分 胁迫 因子 。 蒜 腾 量 与 蒸发 量 较 难 
进行 准确 区 分 ,实际 落 腾 计算 不 准确 将 引起 光合 


表 6 作物 生长 及 产量 模拟 结果 及 评价 


Tab.6 Crop growth and yield simulation results and evaluation 


模型 kı kz Ài CVF Ar r 平均 RE/% 样本 数 
VG 0.6860 0.1143 0.6456 0.5140 0.5185 0.9803 0.9186 8.73 75 
MP 0.8190 0.0725 0.9028 0.6625 0.2374 0.1551 0.9223 8.40 75 
LS 0.8880 0.0409 0.7874 0.6005 0.1140 0.6999 0.9268 8.42 75 


TE CVE 为 光合 产物 转化 效率 ;乒乓 后 人均 为 经 验 系 数 ; 检 验 均 达到 极 显著 水 平 。 


IF @VG ITF (b) MP 9 
$ BE $ 8 $ 8 
之 7 上 è 7 è 7 
i a ql 
BS 6 上 A 6 BS 6 
aK ax a 

5 y=1.003x 5 y=0.9902x 5 y=0.9906x 

e R?=0.7334 R=0.8045 R=0.8151 
A eie u A Le i 1 1 1 Ales 
4 > 6 7 8 9 4 5 6 7 8 4 5 6 7 8 9 


实测 值 /hm 


实测 值 /hm 


图 2 产量 模拟 值 与 实测 值 对 比 


实测 值 /t-hm? 


Fig. 2 Comparison of output simulation value and actual measurement value 
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相 比 作为 影响 因子 ,如 在 FAO 第 33 号 灌 排 文 件 中 给 
出 的 产量 与 水 分 响应 转换 关系 中 采用 了 实际 蒸 散 
量 与 潜在 蒸 散 量 的 比值 ,这 种 经 验 关系 使 用 较为 广 
Schaar 

植物 光合 产物 分 配 系数 模式 ”的 建立 可 实现 
生物 产量 在 根 、 茎 . 叶 籽粒 间 的 分 配 ,从 而 完成 作 
物 生 长 动态 及 产量 的 模拟 。 生 物产 量 在 根茎. 叶 、 
籽粒 间 分 配 时 也 会 受到 水 分 养分 等 环境 因素 的 影 
响 , 即 影响 植株 各 器 官 的 生长 发 育 ,这 种 影响 又 会 
反馈 到 潜在 蒸 散 量 与 潜在 生物 量 模拟 中 ,对 此 张 淑 
杰 等 “建立 了 动态 作物 系数 计算 方法 取代 了 原本 
的 固定 作物 系数 。 本 研究 着 重 对 比 水 分 胁迫 响应 
函数 形式 ,没有 考虑 作物 系数 动态 变化 对 作物 生长 
和 产量 模拟 的 影响 ,该 处 理 对 产量 模拟 没有 造成 过 
大 误差 ,如 康 国芳 汰 模拟 的 冬小麦 产量 平均 相对 误 
差 为 7.10%; 王 伟 等 "9 采用 CERES-Wheat 模型 预测 
冬小麦 产量 ,4 a 平均 相对 误差 为 9.42%, 同 本 研究 
给 出 的 相对 误差 (8.40%~8.73% ) 接 近 。 

不 同 的 水 分 胁迫 响应 函数 形式 没有 对 产量 模 
拟 产 生 明 显 的 影响 ,主要 是 因为 水 分 胁迫 系数 、 养 
分 胁迫 系数 的 参数 反 演 ,在 一 定 程度 上 修正 了 蒸 散 
E . 吸 氮 量 对 作物 生长 动态 及 产量 的 影响 。 利 用 反 
演 模 型 推 求 土壤 水 分 特征 参数 已 成 为 近 些 年 的 研 
究 热点 ”5 ,相对 室内 试验 的 方法 ,参数 反 演 获得 的 
参数 在 一 定 程度 上 概 化 了 不 同 空间 点 的 差异 。 另 
一 方面 ,参数 反 演 对 作物 生长 模型 中 非 保 守 参 数 的 
率 定 也 有 很 大 作用 ,通过 率 定 可 以 使 模型 本 地 化 ， 
显著 提高 模拟 精度 。 


5 结论 


(1) KIA VG MP 及 LS 都 可 以 获得 较 好 的 土壤 
水 分 模拟 结果 及 产量 模拟 结果 ,相关 性 系数 值 均 在 
0.71 以 上 ,土壤 含水 率 模拟 值 与 实测 值 的 平均 RE 小 
F 12.5% ,模拟 产量 的 平均 RE 小 于 8.8%。 

(2) 就 土壤 水 分 动态 模拟 效果 而 言 ,$ 型 水 分 胁 
迫 响应 函数 好 于 四 凸 型 水 分 胁迫 响应 函数 ,VG 水 
分 胁迫 响应 函数 又 好 于 LS 水 分 胁迫 响应 函数 。 采 
用 VG 水 分 胁迫 响应 函数 得 到 的 薰 散 量 更 能 反映 实 
际 情 况 。 


(3) 本 研究 中 提出 的 LS 水 分 胁迫 响应 函数 委 
好 的 避免 了 分 段 函数 中 速率 突变 问题 具有 较 好 的 
模拟 效果 ,但 是 其 随 环境 变化 的 稳定 性 还 有 待 进 一 
步 研究 。 


参考 文献 (References ) 


[1] FSA. 黄土 高 原 林 地 深层 土壤 根系 吸水 过 程 及 其 对 水 分 胁迫 
和 土壤 碳 输 入 的 影响 [D]. 杨凌 : 西北 农林 科技 大 学 , 2019. [Li 


Huijie. Root water uptake process in deep soil for forest growing 


on the Loess Plateau and its effect on water stress and soil carbon 
input[D]. Yangling: Northwest A & F University, 2019. ] 

[2] 康 绍 忠 , 刘晓明 , 能 运 章 . A) AAR RIKARE]. 西北 
农林 科技 大 学 学 报 (自然 科学 版 ), 1992, 20(2): 5-12. [Kang Sha- 


ozhong, Liu Xiaoming, Xiong Yunzhang. Research on the model of 


mater uptake by winter wheat root system[J]. Journal of Northwest 

A & F University (Natural Science Edition), 1992, 20(2): 5-12. | 
[3] ZM. 东北 玉米 物质 生产 及 根系 吸水 对 水 分 胁迫 的 响应 研究 

[D]. 沈阳 : 沈阳 农业 大 学 , 2019. [Jiang Peng. Northeast corn ma- 


terial production and root water absorption study on the response 


of water stress[D]. Shenyang: Shenyang Agricultural University, 

2019. | 

[4 EEH, KEI. HEL He R KERER. 草 业 学 报 ， 
2017, 26(3): 214-225. [Wang Yuyang, Chen Yapeng. Research 


progress in water uptake models by plant roots[J]. Acta Pratacul- 
turae Sinica, 2017, 26(3): 214-225. ] 

[S] EER, 孙 西 欢 , 马 娟 娟 , 等 . 植物 根系 吸水 研究 中. 山西 水 利 ， 
2007(1): 85-86, 88. [Wang Chunxia, Sun Xihuan, Ma Juanjuan, et 


al. Research on plant root water absorption|J]. Shanxi Water Re- 
sources, 2007(1): 85-86, 88. ] 

[6] DRX, 黄 俊 梅 , 王 玉 潮 , 等 . 植物 根系 吸水 特性 研究 由 . 干旱 区 
地 理 , 1999, 22(2): 88-96. [Zhao Chengyi, Huang Junmei, Wang 


Yuchao, et al. Research on water absorption characteristics of 
plant roots[J]. Arid Land Geography, 1999, 22(2): 88-96. ] 

[7] EW. 农业 水 文 模型 中 的 关键 参数 对 作物 蒸腾 量 影响 的 数值 研 
究 [D]. 杭州 : 浙江 大 学 , 2020. [Li Cong. Numerical study on the 


effects of key parameters of agro-hydrological model on crop tran- 
spiration[D]. Hangzhou: Zhejiang University, 2020. | 

[8] AAR AB, FEM, 虽 海 深 . 荒漠 植物 根系 吸水 的 数学 模型 中 . F 
晶 区 资源 与 环境 , 2001, 15(2): 76-80. [Zhu Yonghua, Wu Yanq- 


ing, Lii Haishen. Mathematical model of water absorption of eremo- 
phyte root system|J]. Journal of Arid Land Resources and Environ- 
ment, 2001, 15(2): 76-80. ] 

[9] Feddes R A, Kowalik P J, Malinka K K, et al. Simulation of field 
water uptake by plants using a soil water dependent root extraction 
function|J]. Journal of Hydrology, 1976, 31(1—2): 13-26. 

[10] 金 欣欣 , 石 建 初 , 李 森 , 等 . 根系 吸水 模型 模拟 覆 膜 旱 作 水 稳 气 
FLEE]. 农业 工程 学 报 , 2017, 33(9): 107-115. [Jin Xinxin, 


Shi Jianchu, Li Sen, et al. Modeling stomatal conductance using 


202204.00136v1 


chinaXiv 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


王 铁 英 等 : 根系 水 分 胁迫 响应 函数 对 
root-water-uptake in ground cover rice production system|J]. Trans- 
actions of the Chinese Society of Agricultural Engineering, 2017, 
33(9): 107-115. | 

Bouten W, Heimovaara T J, Tiktak A. Spatial patterns of through- 
fall and soil water dynamics in a Douglas fir stand[J]. Water Re- 
sources Research, 1992, 28(12): 3227-3233. 

Feddes R A, Kowalik P, Zarandy H. Simulation of field water use 
and crop yield[M]. Wageningen: Centre for Agricultural Publish- 
ing and Documentation, 1978: 189. 

吉 喜 斌 , 康 尔 泗 , 陈 仁 升 , 等 . 植物 根系 吸水 模型 研究 进展 由. 西 
北 植物 学 报 , 2006, 26(5): 1079-1086. [Ji Xibin, Kang Ersi, Chen 
Rensheng, et al. Research advances about water-uptake models by 
plant roots[J]. Acta Botanica Boreali-Occidentalia Sinica, 2006, 26 
(5): 1079-1086. | 

吴 训 . LHR SRO EZ SEAR BS EL 2] BO ee AIL 
[D]. 北京 : 中 国 农 业 大 学 , 2019. [Wu Xun. Effects o 


deficit on crop water consumption by transpiration and their quan- 


soil water 


titative delineations[D]. Beijing: China Agricultural University, 
2019. | 

Wu X, Zuo Q, Shi J, et al. Introducing water stress hysteresis to 
the Feddes empirical macroscopic root water uptake model[J]. Ag- 
ricultural Water Management, 2020, 240: 106293, doi: 10.1016/j.ag- 
wat.2020. 106293. 

HENE, 毛 晓 敏 , 雷 志 栋 , 等 . 土壤 水 分 动态 模拟 模型 及 其 应 用 
[M]. 北京 : 科学 出 版 社 , 2009: 65-121. [Shang Songhao, Mao 


Xiaomin, Lei Zhidong, et al. Dynamic simulation model of soil 
moisture and its application[M]. Beijing: Science Press, 2009: 65- 
121. ] 

KYI, 石 建 初 , 左 强 . 基于 作物 水 分 关系 改进 土壤 水 分 胁迫 修正 
系数 的 反 求 方法 四. 水 利 学 报 , 2020, 51(2): 212-222. [Wu Xun, 


Shi Jianchu, Zuo Qiang. Improving the inverse method to estimate 


the soil water stress reduction function based on crop-water rela- 
tions[J]. Journal of Hydraulic Engineering, 2020, 51(2): 212-222. | 
孙 扬 越 , 申 双 和 . (EAE KE EE. 中 国 农业 气 
Z, 2019, 40(7): 444-459. [Sun Yangyue, Shen Shuanghe. Re- 


search progress in application of crop growth models[J]. Chinese 
Journal of Agrometeorology, 2019, 40(7): 444-459. ] 

王 锐 , 李 亚 飞 , 张丽娟 , 等 . 土壤 湿度 驱动 WOFOST 模 型 及 其 适 
MHJ. 中 国 农业 气象 , 2015, 36(3): 263-271. [Wang Rui, Li Ya- 


fei, Zhang Lijuan, et al. WOFOST model based on soil moisture 
driven and its adaptability[J]. Chinese Journal of Agrometeorology, 
2015, 36(3): 263-271. | 

郎 婷 婷 , MBSE Se, 吴 风 华 , 等 . AEF DSSAT REA EY oe ESE DX. E 
要 农作物 用 水 分 析 町 . 干旱 地 区 农业 研究 , 2019, 37(5): 235- 
242, 248. [Lang Tingting, Hao Mengmeng, Wu Fenghua, et al. 


Study on water requirements of major crops in Beijing-Tianjin-He- 
bei region using DSSAT model[J]. Agricultural Research in the Ar- 
id Areas, 2019, 37(5): 235-242, 248. | 

XU DE, SBAI, FIRR. 土地 生产 力 模型 (PS123) 的 校正 和 验证 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


培 水 及 作物 生长 


ChinaXiv 合 作 期 刊 


动态 和 产量 模拟 影响 的 研究 


UL 中 国 农业 气象 , 2005, 26(3): 150-154. [Liu Hong, Guo Wenli, 
Yu Zhenrong. Validation and calibration of land production simu- 
lation model (PS123)[J]. Chinese Journal of Agrometeorology, 2005, 
26(3): 150-154. | 

杨 丽 霞 , 王 仰 仁 . PS123 模型 用 于 山西 省 冬小麦 水 管理 的 适用 
性 评价 中. 灌溉 排水 学 报 , 2014, 33(4/5): 409-413. 
Wang Yangren. Evaluation on the applicability of PS123 model in 


<A 


Yang Lixia, 


Shanxi Province[J]. Journal of Irrigation and Drainage, 2014, 33(4/ 

5): 409-413. ] 

王国 帅 , 史 海 滨 , 李 仙 岳 , 等 . 基于 HYDRUS-1D 模型 的 荒漠 绿 
Wok thie BRAS VERS). 农业 工程 学 报 , 2021, 37(8): 87-98. 


[Wang Guoshuai, Shi Haibin, Li Xianyue, et al. Simulation and 


evaluation of soil water and salt transport in desert oases of Hetao 
Irrigation District using HYDRUS-1D modellJ]. Transactions of the 
Chinese Society of Agricultural Engineering, 2021, 37(8): 87-98. | 
王 世 明 , 范 敬 龙 , 赵 英 , EO TAR AE PS AIT PRE 
土壤 水 盐 运 移 数值 模拟 加 . 干旱 区 地 理 , 2021, 44(4): 1104- 
1113. [Wang Shiming, Fan Jinglong, Zhao Ying, et al. Numerical 


simulation of water and salt migration in desert soil in the lower 
reaches of Tarim River under salt-water irrigationlJ]. Arid Land Ge- 
ography, 2021, 44(4): 1104-1113. ] 

Hansen S, Abrahamsen P, Petersen C T, et al. Daisy: Model use, 
calibration, and validation[J]. Transactions of the ASABE, 2012, 
55(4): 1315-1333. 

胡 克 林 , 深 浩 . 农田 土壤 -作物 系统 过 程 模型 及 应 用 [M]. 北京 : 
科学 出 版 社 , 2019: 1-49. [Hu Kelin, Liang Hao. Farmland soil: 


Crop system process model and application[M]. Beijing: Science 
Press, 2019: 1-49. ] 

姜 志 伟 , 陈 仲 新 , 周 清 波 , 等 . CERES-Wheat 作物 模型 参数 全 局 
敏感 性 分 析 町 . 农业 工程 学 报 , 2011, 27(1): 236-242. [Jiang Zhi- 


wei, Chen Zhongxin, Zhou Qingbo, et al. Global sensitivity analy- 


sis of CERES-Wheat model parameters|J]. Transactions of the Chi- 
nese Society of Agricultural Engineering, 2011, 27(1): 236-242. ] 

陈 小 民 , 崔 红 , 刘 志 铭 , 等 . 作物 模型 在 我 国 玉米 生产 中 的 研究 
与 应 用 轧 . 玉米 科学 , 2018, 26(3): 115-120. [Chen Xiaomin, Cui 


Hong, Liu Zhiming, et al. Research and application on maize pro- 


duction by crop models in ChinalJ]. Journal of Maize Sciences, 
2018, 26(3): 115-120. ] 

EPRI, 张 银 锁 , FIRR. 基于 作物 生长 模型 的 夏 玉米 灌溉 需求 
4y Wr (J). 作物 学 报 , 2003, 29(3): 419-426. [Mao Zhenqiang， 


Zhang Yinsuo, Yu Zhenrong. Water requirement and irrigation 


scenarios of summer maize production aided by crop growth simu- 
lation model[J]. Acta Agronomica Sinica, 2003, 29(3): 419-426. ] 

胡 克 林 , 李 保 国 , 陈 研 , 等 . 作物 生长 与 土壤 水 氮 运 移 联合 模拟 
的 研究 工 一 一 模型 中 . 水 利 学 报 , 2007, 38(7): 779-785. [Hu Ke- 


lin, Li Baoguo, Chen Yan, et al. Coupled simulation of crop growth 
with soil water-heat-nitrogen transport [: Model[J]. Journal of Hy- 
draulic Engineering, 2007, 38(7): 779-785. | 

王 康 . 非 饱 和 土壤 水 流 运动 及 溶质 运 移 [M]. 北京 : 科学 出 版 社 ， 


202204.00136v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


FHEWY 


[33] 


[34] 


[35] 


[36] 


[37] 


[38] 


[39] 


[40] 


[41] 


2010: 1-109. [Wang Kang. Unsaturated soil water flow movement 
and solute transport|M]. Beijing: Science Press, 2010: 1-109. | 

康 绍 忠 . 土壤 -植物 -大 气 连续 体 水 分 传播 理论 及 其 应 用 [M]. 北 
京 : 中 国 水 利水 电 出 版 社 , 1994: 1-121. [Kang Shaozhong. Soil- 


plant-atmosphere continuum water propagation theory and its ap- 
plication[M]. Beijing: China Water Conservancy and Hydropower 
Press, 1994: 1-121. ] 

王 仰 仁 . 考虑 水 分 和 养分 胁迫 的 SPAC 水 热 动态 与 作物 生长 模 
拟 研究 [D]. 杨凌 : 西北 农林 科技 大 学 , 2004. [Wang Yangren. 


Water, heat transfer and crop growth simulation in SPAC with wa- 


ter and nutrient stress[D]. Yangling: Northwest A & F University, 
2004. | 

雷 志 栋 , 杨 诗 秀 , 谢 森 传 . 土壤 水 动力 学 [M]. 北京 : 清华 大 学 出 
版 社 , 1988: 77-130. [Lei Zhidong, Yang Shixiu, Xie Senchuan. 


Soil water hydrodynamics|M]. Beijing: Tsinghua University Press, 
1988: 77-130. ] 

Novak V. Estimation of soil-water extraction patterns by roots|[J]. 
Agricultural Water Management, 1987, 12(4): 271-278. 

FIRR, FER, MEE. 土地 利用 系统 分 析 [M]. 北京 : 中 国 农 
业 科 技 出 版 社 , 1997: 101-162. [Yu Zhenrong, Wang Jianwu, Qiu 


Jianjun. Land use system analysis[M]. Beijing: China Agricultural 
Science and Technology Press, 1997: 101-162. | 

姚 丽 . 限量 供水 条 件 下 精准 灌溉 技术 集成 效益 分 析 [D]. 天 津 : 
天 津 农学 院 , 2020. [Yao Li. Analysis of integrated benefit of pre- 


cision irrigation technology under limited water supply[D]. Tian- 
jin: Tianjin Agricultural University, 2020. ] 

明道 绪 . 田间 试验 与 统计 分 析 [M]. 第 三 版 . 北京 : 科学 出 版 社 ， 
2013: 39-275. [Ming Daoxu. Field experiments and statistical 
analysis[M]. 3" ed. Beijing: Science Press, 2013: 39-275. ] 
TA, 程 煜 , 张 体 彬 , 等 . 利用 HYDRUS-2D 模拟 膜 下 滴灌 玉 
米 农田 深层 土壤 水 分 动态 与 根系 吸水 四 .干旱 地 区 农业 研究 ， 
2021, 39(3): 23-32. [Ding Yuntao, Cheng Yu, Zhang Tibin, et al. 


Modeling of dynamics of deep soil water and root uptake of maize 
with mulched drip irrigations using HY DRUS-2D{J]. Agricultural 
Research in the Arid Areas, 2021, 39(3): 23-32. | 

ABER. 绿洲 农田 土壤 水 分 运 移 规律 及 其 对 作物 生长 的 影响 DD] 
乌鲁木齐 : 新 疆 农业 大 学 , 2004. [Sheng Yu. The law of soil water 


movement and its influence on growth of crop in oasis farmland 


[D]. Urumqi: Xinjiang Agricultural University, 2004. ] 
杨 诗 秀 , MUA, 陆 秀 文 , 等 . 应 用 土壤 物理 -土壤 水 和 温度 应 


[42] 


[43] 


[44] 


[45] 


[46] 


[47] 


[48] 


= 


JIM]. 北京 : 水 利水 电 出 版 社 , 1984: 123-151. [Yang Shixiu, 
Liu Danren, Lu Xiuwen, et al. Applied soil physics: Soil water and 
temperature application[M]. Beijing: China Water Conservancy 
and Hydropower Press, 1984: 123-151. ] 

SEHR, 米 娜 , 明 惠 青 , 等 . WOFOST IRAIZ HC ECE MT EF 
量 模拟 影响 加. 应 用 气象 学 报 , 2021, 32(1): 52-64. [Cai Fu, Mi 


Na, Ming Huiqing, et al. Effects of improving evapotranspiration 


parameterization scheme on WOFOST model performance in simu- 
lating maize drought stress process[J]. Journal of Applied Meteoro- 
logical Science, 2021, 32(1): 52-64. ] 

孙 仕 军 , 张 琳 琳 , 陈 志 君 , 等 . AquaCrop 作物 模型 应 用 研究 进展 
由. 中 国 农业 科学 , 2017, 50(17): 3286-3299. [Sun Shijun, Zhang 


Linlin, Chen Zhijun, et al. Advances in AquaCrop model research 
and application|J]. Scientia Agricultura Sinica, 2017, 50(17): 
3286-3299. | 

张 淑 杰 , 周 广 胜 , FRF. 基于 涡 度 相关 的 春玉 米 逐 日 作物 系数 
KAR BOGUT]. 应 用 气象 学 报 , 2015, 26(6): 695-704. [Zhang 
Shujie, Zhou Guangsheng, Li Rongping. Daliy crop coefficient of 


spring maize using eddy covariance observation and its actual 
evapotranspiration simulation[J]. Journal of Applied Meteorologi- 
cal Science, 2015, 26(6): 695-704. ] 

康 国 芳 . 河北 省 冬小麦 生长 模拟 与 叶绿素 荧光 参数 分 析 [D]. 保 
定 : 河北 农业 大 学 , 2008. [Kang Guofang. Study on winter growth 


simulation and chlorophyll fluorescence parameter analysis in He- 
bei Province[D]. Baoding: Hebei Agricultural University, 2008. ] 

王 伟 , 黄 义 德 , 黄 文江 , 等 . 作物 生长 模型 的 适用 性 评价 及 冬 小 
麦 产量 预测 员 . 农业 工程 学 报 , 2010, 26(3): 233-237. [Wang 
Wei, Huang Yide, Huang Wenjiang, et al. Applicability evaluation 


of CERES- wheat model and yield prediction of winter wheat[J]. 

Transactions of the CSAE, 2010, 26(3): 233-237. ] 

俞 明 涛 , 张 科 锋 . 基于 HYDRUS-2D 软件 的 土壤 水 力 特 征 参 数 
反 演 及 间接 地 下 滴灌 的 土壤 水 分 运动 模拟 串 . 浙江 农业 学 报 ， 

2019, 31(3): 458-468. [Yu Mingtao, Zhang Kefeng. Identification 


of soil hydraulic parameters based on HYDRUS-2D software and 
simulation of soil water movement under indirect subsurface drip irri- 
gation[J]. Acta Agriculturae Zhejiangensis, 2019, 31(3): 458-468. | 

刘 彬 彬 . 基于 微 遗 传 算法 的 土壤 水 力 特征 参数 反 演 及 水 分 运动 
模拟 DD]. 杭州 : 浙江 大 学 , 2018. [Liu Binbin. Identification of soil 


hydraulic property parameters based on the micro- genetic algo- 


rithm and simulation[D]. Hangzhou: Zhejiang University, 2018. ] 


202204.00136v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


王 铁 英 等 : 根系 水 分 胁迫 响应 函数 对 土壤 水 及 作物 生长 动态 和 产量 模拟 影响 的 研究 


Effect of root water stress response function on 
soil water, crop growth dynamics and yield simulation 


WANG Tieying, WANG Yangren', CHAI Junfang, GUO Wenjun 
(1. School of Hydraulic Engineering, Tianjing Agricultural University, Tianjing 300392, China; 


2. Irrigation Experiment Station of Huoquan Water Conservancy Affairs Center, Hongdong County, Linfen 031600, Shanxi, China) 


Abstract: To investigate the effect of the root moisture stress response function (RMSRF) on soil moisture 
dynamic and yield simulation, the soil moisture dynamic, and wheat yield were simulated based on the Richards 
equation and PS123 crop growth model. The three RMSRF of VG (S-curve), MP (convex curve), and LS (S-curve) 
were compared. The soil characteristic and RMSRF parameters in the model were inverted to determine the 
optimal parameter values using the measured data from Huoquan Station (three years) and Xiaohe Station (two 
years) in Shanxi Province, China to obtain the simulation results of soil moisture content, evapotranspiration, and 
grain yield. The analysis is as follows. (1) The soil moisture and yield simulations under VG, MP, and LS 
functions are good. The RMSE value is between 0.021 and 0.036; the r test between the simulated and measured 
values has reached an extremely significant level. For the Huoquan Station, under the three RMSRF of VG, MP, 
and LS conditions, the average relative errors of the simulated and measured values for soil moisture content are 
6.37%, 8.26%, and 7.18%, respectively. The minimum certainty coefficient value is 0.7814. The average relative 
error values of the simulated wheat yield are 8.73%, 8.40%, and 8.42%, respectively. For the Xiaohe Station, the 
average relative error of the simulated and measured values for soil moisture content is greater than that of the 
Huoquan Station, with a maximum of 12.47%. (2) The simulation of soil moisture dynamic shows that RMSRF 
by the S-shaped curve shows better results than RMSRF by the concave-convex curve. The average relative error 
using the VG function is smaller, and the simulated evapotranspiration is closer to the actual measurement. (3) In 
terms of yield simulation results, there is no significant difference among the three RMSRF. In summary, the VG 
function is a high-precision RMSRF, and the model is simple and convenient. The suggested LS function can 
improve the simulation accuracy, but the stability of the model needs further research. 

Key words: moisture stress; root moisture uptake model; farmland moisture simulation; yield simulation; pa- 


rameter inversion 


